Faraday Rotation Spectroscopy of Quantum-Dot Quantum Wells 
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Time-resolved Faraday rotation studies of CdS/CdSe/CdS quantum-dot quantum wells have re- 
cently shown that the Faraday rotation angle exhibits several well-defined resonances as a function 
of probe energy close to the absorption edge. Here, we calculate the Faraday rotation angle from 
the eigenstates of the quantum-dot quantum well obtained with k • p theory. We show that the 
large number of narrow resonances with comparable spectral weight observed in experiment is not 
reproduced by the level scheme of a quantum-dot quantum well with perfect spherical symmetry. 
A simple model for broken spherical symmetry yields results in better qualitative agreement with 
experiment. 

PACS numbers: 78.67.Hc,73.22.-f,78.20.Ls 



I. INTRODUCTION 

Semiconductor heterostructures with a size of several 
nanometers, such as core-shell quantum dots, have been 
widely studied in the past years. The combination of 
several semiconducting materials in a single nanocrystal 
provides additional degrees of freedom which allow one 
to tailor the properties of nanocrystals to a certain ex- 
tent i£ Capping layers of a high-bandgap material on 
a low-bandgap quantum dot (QD) passivate the surface 
and increase the optical gain^4 Quantum-dot quantum 
wells (QDQWs), with several layers of a low-bandgap 
material sandwiched between a high-bandgap core and 
cap are not only of interest from a fundamental point of 
Yw\Y^LL±ii±Llli but could also be scaled up to functional 
heterostructures in a single nanometer-size object. 

Although QDQWs have been studied for several years, 
a detailed investigation of the quantum size spectrum us- 
ing photoluminescence (PL) or absorption spectroscopy 
is challenging because of inhomogeneous broadening. PL 
typically shows a single broad peak, while the absorp- 
tion spectrum exhibits a staircase-like spectrum with- 
out distinct resonances. Individual exciton transitions 
have been observed using techniques such as hole burn- 
ing where a homogeneous subset of QDQWs is probed 
selectively^ Time-resolved Faraday rotation (TRFR), a 
well-established technique to investigate the spin dynam- 
ics in nanocrystalsjii*i£ has recently been used to charac- 
terize colloidal CdS/CdSe/CdS QDQWs with hexagonal 
wurtzite crystal structure J« 3 The TRFR signal amplitude 
was found to depend sensitively on probe energy and 
to exhibit three resonances with linewidths as small as 
10 — 20mcV within an energy window of 0.2 eV around 
the absorption edge. Hence, TRFR does not only provide 
information on the spin dynamics, but also is a sensitive 
spectroscopic technique which allows one to identify in- 
dividual exciton transitions in QDQWs. 

Here, we develop a microscopic theory for the TRFR 
signal amplitude as a function of probe energy for 
QDQWs. From the eigenstates calculated using a two- 
and four-band kp description for the conduction and va- 



lence band states, respectively, we determine the dynamic 
dielectric response functions for circularly polarized 
light and the amplitude of the TRFR signal, 9p(E). 
While several narrow resonances in 0p(E) are predicted 
and the resonance energies are well reproduced by k • p 
theory, our calculations show that the spectral weight 
of the resonances detected experimentally 1 ^ is not repro- 
duced by the level scheme of a spherical QDQW (Sec. lTlJl . 
We discuss a simple model in which deformation of the 
QDQW leads to mixing of valence band multiplets, such 
that spectral weight is re-distributed between the differ- 
ent dipole- allowed exciton transitions. Our calculations 
show that this model yields better qualitative agreement 
with experimental data and may provide a possible ex- 
planation for the well-defined resonances in the TRFR 
signal and the featureless increase in the absorption sig- 
nal (Sec. Illlf) . In Sec. II VI we summarize our results. 



II. SPHERICAL QDQWS 
A. Energy level scheme 

We first consider CdS/CdSe/CdS QDQWs with per- 
fect spherical symmetry as shown in Figs, ^a), (b). 
The CdS/CdSe/CdS QDQWs have hexagonal wurtzite 
crystal structure. For the QDQWs in Ref. HIH the ra- 
dius of the CdS core and the width of the CdS cap 
are r± ~ 1.7 nm and r3 — r-i ~ 1.6 nm, respectively. 
The width of the CdSe quantum well (QW), r 2 — r\ — 
"CdSeacdSo, is estimated from the number of CdSe mono- 
layers, rtcdSc, and the monolayer thickness in bulk CdSe, 
acdSo — 0.43 nm. 14 Because the focus of this paper is the 
TRFR signal amplitude as a function of probe energy, we 
restrict ourselves to the simplest realistic description of 
the QDQW. The energy levels are calculated from k • p 
theory^i^ with a two-band Hamiltonian for the conduc- 
tion band states and a four-band Luttinger Hamiltonian 
for the heavy and light hole valence band states. 

In the two-band approximation, the conduction band 
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FIG. 1: (color online), (a) Schematic representation of the 
CdS/CdSe/CdS QDQW of Ref. [H (b) The radius of the 
central CdS core is denoted by ri, the width of the CdSe QW 
by V2 — r i , and the width of the CdS cap by 7"3 — r% ■ Numerical 
values are given in the text, (c) Radial probability distribu- 
tion of the conduction band ground state 15 e . Because of the 
small conduction band mass, the state is not well localized in 
the QW (indicated by dashed vertical lines), (d) Radial prob- 
ability distribution of the IS3/2 valence band state. The two 
components of the wave function, 7?o and R2, are shown in 
solid and dashed lines, respectively. Because the valence band 
mass is large compared to the conduction band mass, the va- 
lence band states are more strongly localized in the QW. (e) 
Energies of the lS e , 2S e , and lP e conduction band levels as 
a function of QW width, (f) Energies of the lowest three va- 
lence band levels with S^-symmetry (squares) as a function 
of QW width. The energies of the states I-P3/2, I-D5/2, IP5/2, 
and I-D7/2 are also shown (from bottom). Labels for I-D5/2, 
IP5/2, and I-D7/2 are omitted for clarity. 



states are determined by 



2m(r) 



Voir), 



(1) 



where the band mass is given by the value of bulk CdSe 
and CdS, respectively^ 



m(r) 



racdSc = O.llmo, r 1 <r<r 2 , 
mcdS — 0.15 mo, r < r\ orr > r%. 



(2) 



mo denotes the free electron mass and k the momentum 
operator for the envelope wave function. The potential 
V c (r) in i? c b represents the offset of the CdS conduction 
band relative to CdSe, 0.32 eVJ& Vanishing boundary 
conditions are imposed at r$. In Fig. 1(e), the lS e , 2S e , 



and lP e energy levels are displayed as a function of QW 
width. Figure ^c) shows the radial probability distri- 
bution for the conduction band ground state, lS e , for 
^CdSe = 3. Because of the small conduction band mass, 
the state is not well localized in the QW. The lattice 
mismatch at the CdS/CdSe interface (4%) is expected 
to also modify the electronic band structure^ an effect 
that is neglected in our calculations. While the detailed 
band structure close to the interface is not fully under- 
stood, a likely scenario is that the lattice constant varies 
gradually, giving rise to band realignment and a grad- 
ual change in the radial potential V c (r) rather than the 
step-like behavior considered here. 

In the spherical approximation, heavy and light hole 
valence band states are determined by the Luttinger 
Hamiltonian 



H v b = [71 (t 



7(7-) 

to 



(k-J) 2 + K(r) (3) 



for r 7^ rx,2- J denotes the spin operator of the J = 3/2 
valence band multiplet. States of the J = 1/2 split-off 
band are offset in energy by 0.4 eV and are neglected in 
the following. 71 (r) and 7(7") are the Luttinger parame- 
ters of CdSe (CdS) for r inside (outside) the QW| 17 i 20 i 21 



and 



7iW 



7(r) 



7i,CdSc = 1-67, r x <r <r 2 ; 
7i, CdS = 1-09, r<r 1 orr>r 2 , 



7cdSo = 0.56, ri < r < r 2 ; 
7CdS = 0.34, r < r\ orr > r 2 . 



(4) 



(5) 



Because k does not commute with 71 (r) and 7(7"), op- 
erator ordering in Eq. @ is important for r = 7*1,2 (see 
Ref. ^| and references therein) and will be specified by 
the boundary conditions at the interfaces in Eq. JSJ be- 
low. V v (r) represents the offset of 0.42 eV of the CdS 
valence band edge relative to CdSe^ 

The eigenstates of ff v b have been explicitly calculated 
for both QDsS 2 * 2 ^ and QDQWs-i£J£ While the orbital 
angular momentum L of the envelope wave function is 
no longer a good quantum number, the sum of orbital 
and spin angular momentum, F = L + J, commutes with 
the Hamiltonian, [_ff v b,F] = 0. Eigenstates of H v t, are 
labelled according to the quantum numbers F, F Zl and 
the smallest angular momentum component L of the en- 
velope wave function. The eigenstates^ 2 * 2 ^ 



■ 2 , 2 i F, 



\nL F ; F z ) = R L (r)\L, §, F, F z ) + R L +2(r)\L 

(6) 

are superpositions of envelope functions with angular mo- 
mentum L and L + 2. For r < r\, r\ < r < r 2 , and 
r 2 < r, the radial wave functions Rl{t) and Rl+ 2 (t) are 
solutions of the differential equations 2 ^ 



3 



-&M + ^) ( d r + g A - ^) + v v (r) c 275 £ (a 2 + ^ + « j 1 , flL 



C27 



& (a? - ^±ia r + ±e«i) -^h, - c l7 ) (a? + fa, - <i±2<i±a) + v,(r) ) \ 



The full Hamiltonian with the correct operator ordering 
at r = ri 2 is given, e.g., in Ref.ll6land is omitted here for 
brevity. The dimensionless constants Ci )2 depend on F 
and L and can be read off from the values C*i. 2 in Table I 
of Ref.l22l with the correspondence ci 9 = 2C\^/^- Equa- 
tion is solved with a piecewise ansatz in the spherical 
Bessel functions and n^. The boundary conditions for 
the derivatives require that 



71 d r + C17 (d r + ^r) 
-027(^-7) 



-c 27 (a, + ^) 

71 9 r - C17 (<9 r + ^ 

Rh+2 



(8) 



is continuous at ri j2 ^ 

We have calculated the lowest conduction and valence 
band energy levels for the QDQWs of Ref. 0. The 
wave functions of the conduction band ground state, 
XS e , and the components i?o,2 of IS3/2 are shown in 
Figs.^c) and (d), respectively. The lowest valence band 
states are (in order of increasing energy) IP3/2, 1 5*3/2, 
lD 5/2 , and lP 5/2 [Fig. Hff)]. Similarly to CdS/HgS/CdS 
QDQWsji^iS^iiS^ the envelope function of the valence 
band ground state is p-type, implying a dark exciton 
ground state. While the conduction band states are 
strongly delocalized over core, QW, and cap, the valence 
band states are much better localized in the QW because 
of the larger valence band mass. 

The electron-hole Coulomb attraction is calculated in 
first order perturbation theory because the QDQWs are 
in the strong-confinement regime. We neglect the small 
difference in dielectric constants for CdS and CdSe and 
set e = ecdSc = ecds = 9. For an electron in the conduc- 
tion band ground state, lS e , the exciton binding energy 
U is 



U = 



dr e dr h 



' e' h 



47ree J "' " max(r e ,r h ) 

x[R L (r h ) 2 + R L +2(r h ) 2 }, 



<AlS>e) 2 (9) 



with %pis e (r e ) the radial wave function of lS e and Rl,l+2 
the radial components of the hole state. For both the 
I-P3/2 an d the IS3/2 multiplet, the Coulomb integral is 
well approximated by U = — e 2 /47reeo[(Vi + r 2 )/2]. The 
calculation of the dynamic dielectric response functions 
in Sec. Ill Bl involves virtual transitions from one-exciton 
to bi-exciton states, such that the bi-exciton shift must 
be evaluated. For a bi-exciton with two electrons in !S e 



and two holes in 153/ 2 , the bi-exciton shift is of order 
5 meV, sufficiently small that it can be neglected in the 
following. Similarly, the electron-hole exchange interac- 
tion is neglected because its characteristic energy scale is 
of order 1 meVi2& 

The lattice anisotropy of the wurtzite crystal structure 
is taken into account by the anisotropy Hamiltonian 23 



i? an = A (3/2) 2 -J 2 



(10) 



Because the hole wave functions are localized predomi- 
nantly in the CdSe QW, we approximate A = 25meV 
by the CdSe bulk valuei2i2£ H an lifts the degeneracy of 
an Li?-multiplet and splits it into (2F+ l)/2 energy dou- 
blets. In particular, for S3/2 multiplets, the energy shifts 
induced by H an are2£ 

A n>3 /2 = (7iS3/2;3/2|£ an |7iS 3 /2;3/2) ( lla ) 

= A| / drr 2 Rl 



A n>1/2 = <7iS3 /2 ;l/2|if an |nS 3/2 ;l/2) 

1 



(lib) 



drr 2 Rl 



dr r 2 Rn 



for \F Z \ = 3/2 and \F Z \ = 1/2, respectively. 



B. Faraday rotation angle 

We next calculate the amplitude of the Faraday ro- 
tation (FR) angle as a function of probe energy, 0p(E), 
taking into account the single-particle levels, Coulomb in- 
teraction, and crystal anisotropy as described in Sec. Ill Al 
In the pump-probe scheme of Ref. Il 'ii a er~ circularly po- 
larized pump pulse with energy large compared to the ab- 
sorption edge excites spin-polarized excitons. While the 
conduction band electron typically retains its spin po- 
larization on relaxation to the conduction band ground 
state, lS^, the hole spin is believed to randomize quickly 
during relaxation to l-P3/ 2 - The net spin polarization 
along the pump direction, which results from conduction 
band electrons in the spin state \s z =|), is experimen- 
tally detected with TRFR. For probe energy E, 9 F (E) 
is proportional to the difference of the dynamic dielec- 
tric response functions for circularly polarized light, 
e+(E/H) — €-(E/H), which are determined by the optical 



4 



dipole transition matrix elementSf22i22^2i2i 

e F {E) = CE v\(XX\p x + aip y \X in )i 

ar=±l;\XX) 

E — (Exx — Ex in ) 



XX 



Ex ia f+7xx 



(12) 



\Xia) denotes the initial lS e — IP3/2 exciton prepared by 
the pump pulse with an \s z =|)-electron in lS e . The sum 
extends over all bi-exciton states \XX) with energy Exx, 
and "/xx denotes the linewidth of the corresponding bi- 
exciton transition. For simplicity, we assume that the 
linewidth of all transitions is equal, jxx — r. The con- 
stant C depends on the sample size and refractive index. 
6p(E) is finite because transitions to bi-exciton states 
with both conduction band electrons in lS e are allowed 
only if the electrons form a singlet state. Hence, the ma- 
trix elements in Eq. i|12|l can be expressed in terms of the 
transitions from an arbitrary valence band state |$„) to 
the unoccupied conduction band state |lS e ; |), 



9 F (E) = CE <r\(lS e ;l\p x + aip y \$ v )\ 

c7=±l;|*„) 

E — Ex,v 



X ■ 



(E-E x , v ) 2 +T 2 



(13) 



with Ex,v the energy of the lS e -<S> v exciton. As men- 
tioned in Sec. Ill Al the bi-exciton shift is negligible. This 
expression depends only on the single-particle levels in 
the conduction and valence band and can be evaluated 
from the wave functions calculated above. 

The dipole transition matrix elements are evaluated 
following Ref. The overlap integral of the enve- 

lope wave functions is finite only for s-type valence band 
states, |$ ) = \nS 3 / 2 ;F z ). We first consider a QDQW 
with crystal symmetry axis aligned along the laser direc- 
tion. Then, 



\{lS e ;l \p x - ip y \nS 3 /2] F z )\* 
_ 2 
~~ 3 



(14a) 



drr 2 i) ls MS\p z \Z) 



- ^(m Q V) 2 I drr 2 i; 1Sc R 
\{lS e ;l \px + ip y \nS 3/2 ;F z )\ 2 

drr 2 iP 1Ss R (S\p z \Z) 

2(m Q V) 2 [ drrVis e -Ro 



(14b) 



-3/2 



-3/2i 



while the matrix elements vanish for F z = 3/2,-1/2. 
The Kane interband matrix element (S\p z \Z) varies spa- 
tially for the QDQW, but because the valence band 
states are well localized in the CdSe QW, we approxi- 
mate |(S'|p z |Z)| ~ rrioV by the CdSe interband matrix 
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TABLE I: Radial overlap integrals, Eq. H15II . for IS3/2, 2S 3 / 2 , 
and 3S3/2 valence band states as a function of the QW width, 
ncdSe- 
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FIG. 2: (color online), (a) Schematic representation of the 
transitions which contribute to the FR angle, 9f(E), in the 
spectral representation Eq. 1131 . The dipole transition ma- 
trix elements with lS e are finite only for 5'3/2-multiplets. 
(b) 9f(E) calculated from Eq. 1161 for ncdSc = 3 and 
r = 15meV. nS-i/2 multiplets with n > 4 have been ne- 
glected. The transitions from (a) are indicated by arrows. 



element. The overlap integral in Eq. (|14f> . in the follow- 
ing denoted by 



In 



drr 4>is c Ro 



(15) 



depends on the quantum number n via R . For IS3/2, 
2S3/2, and 3S3/2, the numerical values are given in Ta- 
ble [I] for different QW widths. 

With the transition matrix elements, 0p(E) is readily 
evaluated for E close to the absorption edge. We define 
the energy E x , n = Eis e + E nSa/2 +E g + U of the 15 e - 
nS 3 /2 exciton, where E g ~ 1.75 eV is the bandgap of bulk 
CdSe at room temperature and U the Coulomb integral 
[Eq. ©]. Combining Eqs. JT3) and (d, 



C 



E - {E X ,n + A„ !3 / 2 ) 



[E - (Ex, 
E — (Ex, n ■ 



+ A„, 3 / 2 )] 2 
A n ,i/ 2 ) 



3 [E - (E x ,n + A„, 1/2 )] 2 



r 2 



(16) 



where C = 2C(moV) 2 and E was assumed to be suf- 
ficiently close to the absorption edge that transitions 
from the split-off band can be neglected. The energy 
shifts A„ j3 / 2 and A^i/2 are induced by the structural 
anisotropy [Eq. I|ll|)]. 

For every valence band multiplet nS 3 /2, Of{E) shows 
a double-resonance structure, with a main resonance cor- 
responding to transitions from \F Z \ = 3/2 [first term 
in the square bracket in Eq. Ijlbj) ] and a second reso- 
nance with 1/3 smaller spectral weight that corresponds 



5 



to transitions from \F Z 



1/2 states and is shifted 



by A„ x /2 — A„ 3 /2- The second resonance can only 
be resolved if T < A. For the colloidal QDQWs in 
Ref. [H the inhomogeneous line broadening is substan- 
tial, r > 15 meV, and the splitting of resonances caused 
by crystal anisotropy cannot be resolved. The lowest 
transitions with finite matrix elements in the spectral 
representation in Eq. I|13[l are shown schematically in 
Fig-IHa). Figure |2Ib) shows 9f{E) close to the absorp- 
tion edge for ncdSc = 3. While resonances corresponding 
to IS3/2, 25*3/2, and 35 3 / 2 can be resolved, the spectral 
weight for 15*3/2 is significantly larger than that for 25 3 /2 
and 35 3 / 2 - The behavior in Fig. |2Jb) is clearly distinct 
from experimental datai^ which exhibit three resonances 
with comparable spectral weight within 0.2 eV of the ab- 
sorption edge. 

We next calculate 9p (E) for the case that the QDQW 
symmetry axis is not aligned with the laser. By averag- 
ing over different relative orientations, we also determine 
9p{E) for an ensemble of randomly oriented QDQWs. If 
the QDQW crystal symmetry axis is tilted by ifi relative 
to the laser propagation direction, the electric field com- 
ponents of the pump and probe laser pulses, given in the 
laboratory coordinate system (indicated with a super- 
script L in the following), are rotated relative to the spin 
quantization axis of the QDQW. The spectral represen- 
tation for 9p(E) is still given by Eq. (|13|) . but the matrix 
elements must be modified to (15 e ; | L \p x + aipy\<& v ) to 
account for the fact that the pump and probe pulses are 
defined in the laboratory coordinate system L. To eval- 
uate the matrix elements, operators in the L-frame are 
expressed in terms of the QDQW coordinate system, | [ L 
) = isin(V>/2)| T)+cos(>/2)| 1), p£ = cos(ip)p x +sin(ip)p z , 
and py = p y . With the convention for the valence band 
basis functions in Refs. 



(15 e ;4 L \p L x ±ip^\nS z/2 -F z ) = y^j2m V 



x < 



-isin(V>/2)[cos(V>) =F 1], 

-(2/a/3) sinO/2)sin(?A) 
f(i/V3) cos(W2)[cos(V>)Tl], 

(2/V3) cos(W2) sin(V>) 
-(i/\/3)sin(V72)[cos(V)±l], 

icosO/2)[cos(-0) ± 1], 



F z = 
F z = 

F z = 

F z = 



(17) 
3/2; 
1/2; 

-1/2; 
-3/2. 



Inserting these transition matrix elements into the 
spectral representation in Eq. I|13l) . we find that 
the expression for 6p(E) only acquires a pre- factor 

cos 2 (^/2) cos(V'), 



9 F (E) ~ C"^cos 2 (W2) cos(V>) 



(18) 



E- (E X .n + A n , 3/2 ) 
[E - (E x , n + A n , 3 / 2 )] 2 + r 2 

1 E — (Ex,n + A„ (1 / 2 ) 



i[E-{E x ^ + ^ l/2 )f+T^ 



but the relative weight of the individual terms is not al- 
tered. In particular, Eq. l(T%|) implies 9p = for ijj = ir/2. 
This can be understood in physical terms because TRFR 
probes the spin polarization along the QDQW anisotropy 
axis, while, for %j) — 7r/2, the pump pulse only generates 
spin polarization perpendicular to the anisotropy axis. 

Averaging Eq. Ijl8(l over an ensemble of randomly ori- 
ented QDQWs is performed by integration over tj). For 
random QDQW orientation, Eq. (|16f) remains valid with 
the substitution C -> 5C/12. While random QDQW 
orientation reduces the total amplitude of 9p(E), the 
relative spectral weight of the individual contributions 
is not altered and 9p{E) is still as shown in Fig. 0{b). 



III. BROKEN SPHERICAL SYMMETRY 

The calculations in Sec. [H] above show that, for a 
spherical QDQW, 9p(E) exhibits a pair of resonances 
for every 71S3/2 — 15 e exciton transition. However, the 
spectral weight of the resonances decreases rapidly with 
increasing n (Table , such that the large number of 
resonances with comparable spectral weight observed ex- 
perimentally^ is not correctly reproduced by Eq. 11611 . 
The experimental data imply that spectral weight is re- 



distributed from the 15 



3/2 



15 P exciton line to other 



transitions. 

We show next that broken spherical symmetry is a pos- 
sible mechanism which accounts for the experimental FR 
data by mixing of the 153/2 and li"3/ 2 valence band mul- 
tiplets. The admixture of s-type to p-type multiplets re- 
distributes the spectral weight and increases the number 
of resonances with comparable amplitude in 9p(E). As 
will be shown below, the redistribution of spectral weight 
also explains the absence of pronounced resonances in 
the absorption signal. Compared to spherical QDQWs, 
broken symmetry gives rise to a larger energy splitting 
between the lowest valence band states with dominant p- 
type and s-type envelope wave functions, consistent with 
the large Stokes shift between the PL peak and the ab- 
sorption edge observed for CdS/CdSe/CdS QDQWs-i2J£ 

On a microscopic level, broken spherical symmetry 
would most probably result from a spatial variation of 
the QW width [shown schematically in Fig. Eta)]. For 
QWs with ncdSc = 2,3,4,5, a monolayer variation in 
QW width translates into a perturbation with a typical 
energy scale of order 0.1 eV for the hole states. Conduc- 
tion band states are less strongly affected because they 
are not localized in the QW and the energy level spacing 
is much larger. We do not attempt to describe breaking 
of the spherical symmetry microscopically, but restrict 
our discussion to a simple model in which the valence 
band Hamiltonian includes a potential 



5V(r) = v sin 0(1 



(19) 



which mixes 15 3 / 2 and 1^3/2- 9 and (f> denote the az- 
imuthal and polar angle of r relative to the lattice sym- 
metry axis, respectively, vq is a fit parameter. 
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FIG. 3: (color online), (a) Cross section of a QDQW with 
broken spherical symmetry. A deformation of the central CdS 
QD implies a spatial variation of the QW width, which is 
modelled by Eq. 119H . (b) Schematic representation of the 
level scheme for a spherical QDQW and a QDQW with broken 
spherical symmetry. Broken spherical symmetry mixes the 
IS3/2 and IP3/2 multiplets and renders transitions to lS e 
bright for all eight states of the multiplet. 



The exact hole states with the perturbation Eq. i|19|) 
can, at least in principle, be obtained by diagonal- 
ization of 5V(r) in the basis \nLp;F z ) of the spher- 
ical system. The only finite matrix elements are 
(n'L' F ,; F z \SV(r)\nL F ; F z ) oc 8f^fi for \L— L'\ an integer 
multiple of 2 and (n'L'p,; F^\SV(r)\nL F ; F z ) oc 5f„,f^±i 
for \L — L'\ = 1 or \L — L'\ = 3. In order to keep the sub- 
sequent discussion transparent, we restrict the diagonal- 
ization of SV(r) to the subspace spanned by the valence 
band multiplets 153/21 25 3 / 2 , 35*3/2, and 1^3/2- These 
states are dominant for the optical response close to the 
absorption edge and we expect this procedure to qualita- 
tively capture the essential features of broken spherical 
symmetry^ 3 . 

The off-diagonal matrix elements of SV(r) between the 
IS3/2 and IP3/2 multiplet are calculated by expansion of 
the envelope wave functions in spherical harmonics. We 
find 



(1S 3/2 ;F' z \5V(t)\1P 3/2 ;F z ) 



(20) 



x I dri 



4 3 

RqRi + -R2R1 + 7P2P3 

5 5 



for (F z , F z ) = ±(3/2,1/2) or ±(1/2,3/2), and 



(lS 3/2 ;F^\6V(r)\lP 3/2 ;F z 



vq 
v/5 



(21) 



x drr 



2 8 2 
-RqRi + —R2R1 + -R2R3 

3 15 5 



for (F Z ,F Z ) = ±(1/2,-1/2). The diagonal ma- 
trix elements (lS 3 /2] F z \5V(r)\lS 3 / 2 ; F z ) — nvo/i and 
(1P 3/2 ;F z \6V(t)\1P 3/2 ;F z ) ~ nv {9 + 2*| F ,|, 1/2 )/40 are 
evaluated analogously. 

The energy eigenstates are calculated by diagonaliza- 
tion ff v b + H an + SV(r) in the eight-dimensional space 
spanned by the IS3/2 and IP3/2 multiplets. Mixing of 
253/2 and 353/2 with IP3/2 is neglected because of the 
large energy difference [Fig. ^f)]. Because i? v b, Pan, 



and 5V(r) are even under reflection at the x-y-plane, the 
eigensystem consists of four doublets. The eigenstates 



I* 



{SP)s 



3/2 

E 

F,=-3/2 



ai,F z [LS 3 / 2 \F z 



■Pi,FjiP 3 /2;F z )) 

(22) 

have components with s- and p-type envelope functions 
with expansion coefficients ol^f* and fii,F z , respectively. 
The corresponding energy eigenvalue is denoted by E v ^. 
The mixing of 153/ 2 and IP3/2 multiplets renders tran- 

sitions from all eight states i ) to 15 e dipole-allowed 
[shown schematically in Fig. E(b)]. Defining the energy 
of the lS e -¥ v S p exciton by E x% , for QDQWs with sym- 
metry axes oriented along the laser beam, the spectral 
representation for 0p(E) reads 



,(P) ~ C'EhJ2 [\ot 



E 



-3/2 1 
■Ex. 



I«i,l/2| ; 



(23) 



+ C'E In[ 



(E-Ex^+T 2 ' 

E - {E x ,n + A n , 3 / 2 ) 



n=2,3 

1 E 



[E - {E x , n 
(Ex,n + A ?i 



A„,3/ 2 )] 2 + r 2 



,1/2) 



3[p-(p^„ + A n!l / 2 )] 2 + r 2 



( (SP) ) 



For 

the 253/2 and 35 3 /2 valence band multiplets [third 
and fourth line of Eq. l(25ft|. A„ )Fa = A njJ ? z + 
(n5 3 / 2 ; F z \SV(r)\nS 3 / 2 ; F z ) takes into account the energy 
shifts effected by the lattice anisotropy P an and the shape 
anisotropy 5V(r) to first order. For vo = 0, this ex- 
pression correctly simplifies to the one obtained for the 
spherical system. 

In experiment, the QDQWs are randomly oriented rel- 
ative to the laser beams. For QDQWs with symmetry 
axes rotated by ip relative to the laser beam, 9f(E) is 
calculated similarly to Sec. Ill Bl from the matrix elements 
listed in Eq. (|17fl . The experimentally more relevant ex- 
pression for an ensemble of randomly oriented QDQWs is 
found by averaging over ip. For E close to the absorption 
edge, the ensemble average 9f(E) is 



C'Eh 



36 

2 



El 



15 1 a,- 



-3/2 1 



-1/2 1 



(24) 



5|o;i,i/2| 2 - 3|o i:3 / 2 | 2 + 4v3(a»,3/2aj ) i/ 2 



-3/2 c 



,-1/2)] 



E-E 



A', 



(E-Ex,) 2 +T^ 



where o^f* are the expansion coefficients defined in 
Eq. J22J. Transitions involving 25 3 / 2 and 35 3 / 2 are mod- 
ified in an analogous way, but the expressions are omitted 
here for brevity. Note that the relative spectral weight 
of the individual terms changes compared to an oriented 
sample [Eq. ffgjl]. 
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We next discuss how the TRFR signal amplitude, 
Op(E), changes with increasing asymmetry potential ^o- 
FigureOJa) shows the experimental data for 9f(E) (sym- 
bols) and the absorption spectrum (dashed line) for the 
sample with ncdSe = Figures ^b)-(d) show the 
calculated TRFR signal amplitudes (solid lines) for (b) 
t'o = OmeV, (c) v — 40meV, and (d) v — 70meV, 
respectively. The TRFR signal for random QDQW ori- 
entation is calculated from Eq. (|24|l taking into account 
transitions from the 15*3/2, 25 3 / 2 , 35*3/2, and IP3/2 va- 
lence band multiplets. The level broadening T = 15meV 
was chosen to be comparable to the smallest linewidth 
in the experimental data. For the system with spherical 
symmetry, vq — meV, the spectral weight of the reso- 
nance with the smallest energy is significantly larger than 
that of the two resonances at higher energies, in contrast 
to the experimental data. With increasing vq, the spec- 
tral weight is redistributed from the 15 3 / 2 — 15 e tran- 
sition to resonances involving other valence band mul- 
tiplets, such that the spectral weight of the resonance 
with lowest energy approaches that of the higher-energy 
resonances [Figs.^Jc) and (d)]. The corresponding theo- 
retical curves are in better qualitative agreement with the 
experimental results than the TRFR signal calculated for 
a spherically symmetric QDQW. 

The single-particle level spectrum and transition ma- 
trix elements calculated above also allow one to calculate 
the absorption spectrum, 



abs(P) <x 



E 

l*c>;|*„> 



K$c|Ps + ip v |$f 



(25) 

The sum extends over all conduction and valence band 

states and Ex denotes the energy of the corresponding 

exciton transition. Close to the absorption edge, transi- 
ts p) 

tions from |^ i ') and the 25 3 / 2 and 35 3 / 2 multiplets to 
15 e and lP e , respectively, are dominant. The larger num- 
ber of optically allowed transitions in the system with 
broken spherical symmetry [Fig. EJb)] may also explain 
the absence of distinct peaks in the absorption spectrum, 
because transitions from p-type valence band states to 
lP e spectrally overlap with transitions from s-type va- 
lence band states to 15 e . In Figs. Efb)-(d), we show 
the evolution of the absorption spectrum with increasing 
asymmetry of the QDQW (dashed lines) calculated us- 
ing the same parameters as for the TRFR signal. While 
the TRFR signal is dominated by transitions involving 
n5 3 / 2 multiplets fSec. Ill B|) . the absorption spectrum also 
involves transitions from excited valence band multiplets 
such as IP1/2 and 2P 3 / 2 which are neglected here. For the 
absorption spectrum, our calculations restricted to the 
15 3 / 2 , 25 3 / 2 , 35 3 / 2 , and lP 3 /2 multiplets yield valid re- 
sults only within ~ 0.15 cV of the absorption edge. For a 
QDQW with spherical symmetry, pronounced resonances 
are predicted to appear in the absorption spectrum close 
to the absorption edge [Fig. E^b)]. For increasing vq, the 
redistribution of spectral weight effected by mixing of dif- 
ferent valence band multiplets leads to a broadening of 
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FIG. 4: (color online), (a) Experimental data for the ampli- 
tude of the TRFR signal, 6f(E), measured at T = 294 K, and 
the optical absorption as a function of probe energy, (b)-(d) 
Calculated amplitude of the TRFR signal, 9f{E), (solid line) 
and optical absorption (dashed line) for different strength of 
an anisotropy potential, vo = 0, 40, 70meV, respectively. The 
calculation is restricted to the I-P3/2, IS3/2, 2S*3/2, and 35*3/2 
valence band multiplets which dominate the optical response 
close to the absorption edge. For A < 560 nm, transitions 
from higher valence band multiplets become important for 
the absorption spectrum and the restriction to only four va- 
lence band multiplets is no longer valid. Because T = 294 K 
is small compared to the conduction band level splitting, all 
electron spins are assumed to occupy lS e in the calculation 
of F (E). 



these resonances [Fig. life)] which, ultimately, evolve into 
a step-like feature comparable to the experimental data 
[Fig. ffid)]. 

However, for the small transition linewidth T = 
15meV chosen in Fig. Ud), some resonance features can 
still be resolved in the theoretical absorption spectrum. 
One possible explanation for the discrepancy with exper- 
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imental data is that the linewidth of exciton transitions 
to lP e is larger than for lS e because of orbital relaxation 
from lP e to lS e . We also point out that, for vq ~ 70 meV, 
the restriction of our analysis to only four valence band 
multiplets (IS3/2) 2S3/2, 3S 3 / 2 , and I-P3/2) can no longer 
be justified rigorously. A calculation taking into account 
all multiplets in Fig. ^f ) would be required to yield rig- 
orous results for the absorption spectrum in this regime 
of strong symmetry breaking. 

From our analysis of a simple model for broken spher- 
ical symmetry, we conclude that mixing of different va- 
lence band multiplets redistributes the spectral weights 
in the Lehmann representation of 9f (E) and accounts for 
a large number of resonances with comparable spectral 
weight (Fig. E}. Broken spherical symmetry is also con- 
sistent with a large Stokes shift and a featureless increase 
of the absorption spectrum^ 

IV. SUMMARY AND DISCUSSION 

Motivated by recent experiments^ we have cal- 
culated the conduction and valence band states of 
CdS/CdSe/CdS QDQWs using k • p theory. The single- 
particle spectrum allowed us to evaluate the amplitude 
of the TRFR signal as a function of the probe energy, 
9p(E), for samples with well defined and random ori- 
entation of symmetry axes, respectively. For spherical 
QDQWs, 0p(E) exhibits a pair of resonances for every 
nS 3 /2 valence band multiplet [Fig. |2£b)], but does not 
correctly reproduce the spacing and spectral weight of 
the experimentally observed resonances. By contrast, a 
simple model with broken spherical symmetry reproduces 
the experimental data at least qualitatively. 



For both the TRFR signal and the absorption spec- 
trum, the model with broken spherical symmetry yields 
better agreement with experimental data than the cal- 
culations for a spherical QDQW (Fig. 0}. However, the 
theoretical curve does not reproduce all experimental fea- 
tures and large values of vq are required in order to ex- 
plain the experimental data. We, hence, re-emphasize 
that the discussion in Sec. lIIIl is only a qualitative analysis 
of TRFR for broken spherical symmetry. Additional the- 
oretical work is required to improve the quantitative un- 
derstanding of the energy level scheme of CdS / CdSe /CdS 
QDQWs and of the FR angle as a function of probe en- 
ergy 

In particular, more realistic microscopic models for 
the variation of the QW width should be considered. 
Eight-band k • p calculations would allow one to quan- 
tify whether mixing of light and heavy hole states with 
the conduction and split-off band also reduces the spec- 
tral weight of the 153/2 — lSe exciton transition. If the 
full eight-band model gives results comparable to our de- 
scription, this would further corroborate that spherical 
symmetry breaking must be taken into account to under- 
stand the level spectrum of CdS/CdSe/CdS QDQWs. 16 
Tight-binding calculations of the energy spectrum 25 " 34,35 
and 9f(E) would allow one to explicitly include the atom- 
istic structure of interfaces and to compare our calcula- 
tions with a different theoretical framework. 
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